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Abstract 

We present lattice results for baryon number, strangeness and electric charge fluctuations 
as well as their correlations at finite temperature and vanishing chemical potentials, i.e. under 
conditions relevant for RHIC and LHC. We find that the fluctuations change rapidly at the transi- 
tion temperature T c and approach the ideal quark gas limit already at approximately 1.5r c . This 
indicates that quarks are the relevant degrees of freedom that carry the quantum numbers of con- 
served charges at T > l.5T c . At low temperature, qualitative features of the lattice results are 
well described by a hadron resonance gas model. 



1. Introduction 

Fluctuations of conserved charges, like baryon number, electric charge and strangeness, are 
generally considered to be sensitive indicators for the structure of a thermal medium produced in 
heavy ion collisions [1]. Under conditions met in current experiments at RHIC as well as in the 
upcoming heavy ion experiments at LHC the net baryon number is small and QCD at vanishing 
chemical potential provides a good approximation. In this region the transition from the low tem- 
perature hadronic to the high temperature plasma regime is continuous and fluctuations are not 
expected to lead to any singular behavior. Indications for the existence of critical points can only 
show up in higher order derivatives of the QCD partition function with respect to temperature 
or chemical potentials [2]. Through the analysis of fluctuations of conserved charges as well as 
their higher moments and correlations we thus gain insight into the relevant degrees of freedom 
of the system under consideration and at the same time gather information on possible nearby 
singularities in the QCD phase diagram. 

At vanishing baryon number (B), electric charge (Q) and strangeness (S) fluctuations of 
these quantities can be obtained by starting from the QCD partition function with non-zero light 
and strange quark chemical potentials, fL U) d, s = fi u ,d,s/T, The quark chemical potential can be 
expressed in terms of chemical potentials for baryon number (mb), strangeness (fis) and electric 
charge (fi Q ), 
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Hu = jMb + -MQ , Md = - -HQ , Ms = -Mb - -HQ ~ Ms • (1) 

Moments of charge fluctuations, 5Nx = Nx - (Nx), with X = B, Q or S and their correlations are 
then obtained from derivatives of the logarithm of the QCD partition function, i.e. the pressure, 
evaluated at Mb,q,s = 0, 

x b,q,s = d i+ i +k P IT 
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Figure 1 : The quadratic fluctuations of light (left) and strange (middle) quark number versus temperature on lattices of 
size 24 3 X 6 and 32 3 X 8. The horizontal solid lines show results for the ideal massless quark gas; the vertical lines mark 
the temperature interval 185 MeV < T < 195 MeV. The ratio of s and u quark fluctuations (right) are calculated on 
lattices with temporal extent N T = 8. Curves in the right hand figure show results for a hadron resonance gas including 
resonance up to m max = 1.5GeV (upper branch) and 2.5GeV (lower branch), respectively. 



with fix = l-ix/T. While the first derivatives, i.e. baryon number, electric charge and strangeness 
densities, vanish for B,Q,S = 0, their moments and correlation functions with i + j + k even 
are non-zero. The basic quantities we will analyze here are the quadratic, quartic and 6th order 
cumulant of fluctuations, 

*f = ^3<^>' X* = ^fJ«N x )-3(N 2 x ) 2 ), xl = ^3«iVl)-15(^><iVl) + 30<^>3). ( 3 ) 

The gauge field configurations that have been used to evaluate the above observables, had 
been generated previously in calculations of the QCD equation of state jH 01 and the transition 
temperature H,!^. An improved staggered fermion action (p4-action) [7] that strongly reduces 
lattice cut-off effects in bulk thermodynamics at high temperature has been used. The strange 
quark mass has been tuned close to its physical value and the light quark masses have been 
chosen to be one tenth of the strange quark mass. This corresponds to a line of constant physics 
on which the kaon mass is close to its physical value and the lightest pseudo-scalar mass is about 
220 MeV. 

These calculations have been performed on 16 3 x 4, 24 3 x 6 and 32 3 x 8 lattices. This allows 
us to judge the magnitude of systematic effects arising from discretization errors in our improved 
action calculations. The spatial volume has been chosen to be V l ^T = 4, which insures that 
finite volume effects are small. In this calculations, the above observables have been evaluated 
in the temperature interval 0.8 < T/T c < 2.5. 



2. Results 

Before entering a discussion of fluctuations of B, Q and S it is instructive to look into fluc- 
tuations of the partonic degrees of freedom, the light and strange quarks. In Fig. Q]we show 
the temperature dependence of quadratic fluctuations for light quarks^ an d strange quarks;^- 
They change rapidly in the transition region as the carriers of the quantum numbers are heavy 
hadrons at low temperatures but much lighter quarks at high temperatures. In the continuum and 
high temperature limit, these fluctuations quickly approach the Stephan-Boltzmann limit, i.e. an 
ideal massless quark gas. 

The lattice cut-off effects for quadratic fluctuations are shown to be small by comparing the 
measurements on N T = 6 and 8 lattices. Moreover, we have used two different cut-off schemes, 
p4 and asqtad [8]. Both schemes are improved up to O(cr) for thermodynamic quantities at 
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Figure 2: Quadratic, quartic and 6th order cumulant of fluctuations for baryon number, electric charge and strangeness. 
All quantities have been normalized to the corresponding free quark gas values. 



high temperature. Indeed, the quark fluctuations have small deviations from SB limit at high 
temperature on N T = 8 lattices for both schemes. The measurements using p4 and asqtad actions 
are generally in good agreement, although we observe differences in the temperature interval 
between T c and 1.5 T c . 

We observe that the strange quark fluctuations rise slower than the light quark fluctuations. 
As temperature decreases, their ratio x\ix\ (Fig- 03 drops quickly in the transition region. The 
hadron resonance gas (HRG) model calculations show a similar trend of exponential fall at low 
temperature. This can be understood in the zero temperature limit: the light quark fluctuations 
are sensitive to pions, X2IT 2 ~ ex P( -ff Wr), while the strange quark fluctuations are sensitive to 
the lightest hadronic state that carry strangeness, ^/T 2 ~ exp(— m^-/!'). The HRG models can 
not reproduce the lattice measurements at higher temperature, as it breaks down at T 150 MeV, 
where HRG models with different spectrum cuts start to show deviations. 

In Fig. [2]we show results for quadratic fluctuations^, quartic fluctuations^ on N T = 4 and 
6 lattices and 6th order cumulants of fluctuations on N T = 4 lattices, where X — B, Q and S . 
As can be seen, in all cases the quadratic fluctuations rise rapidly in the transition region where 
the quartic fluctuations show a maximum and 6th order cumulants change sign. This generic 
form is indeed expected. In the chiral limit the 4th order cumulants will have a peak and the 6th 
order cumulants will diverge, e.g. the Baryon number cumulants are expected to scale like 

2—n—a 

(4) 

where a —0.25[— 0.015] denoting the critical exponent of the specific heat in 3-d, (9(4) [(9(2)] 
spin models ^ . 

At temperature T > 1.57^ the quadratic, quartic and 6th order cumulants of fluctuations of 
B, Q and S are well described by those of an ideal, massless quark gas. At low temperature, 
we compare lattice results for the ratios of the cumulants with HRG model calculations (Fig. 
[3]). In the framework of a HRG model the ratio of fourth and second order cumulants of baryon 
number fluctuations is independent of the actual value of hadron masses; (xf/Xo)HRG = 1 if all 
hadrons are heavy on the scale of the temperature. This is well reproduced by the lattice results. 
However, in the chiral limit it is expected that the cusp in^* (Fig. [2]) is expected to become more 
pronounced and thus more prominent also in the ration* 

F° r x\lx\ ar, d X4 IXzi tne lattice results also agree with HRG models qualitatively. The 
cumulants of electric charge fluctuations receive contribution from the lightest hadrons, pions. 
In our calculations the masses of the lightest pseudoscalar states, 220 MeV, are close to but not 
exact as the physical pion masses. Therefore, we have analyzed the cumulant ratio X4/X2 m 
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Figure 3: The ratio of fourth and second order cumulants of baryon number (first), strangeness (second) and electric 
charge (third) fluctuations and the ratio of sixth and second order cumulants of electric charge fluctuations (fourth). In 
the last two cases we show curves for a HRG model calculated with physical pion masses (upper curve), pions of mass 
220 MeV (middle) and infinitely heavy pions (lower curve). 

the HRG models with the pion masses at physical value 140 MeV, 220 MeV and infinity. We 
observe that the cumulant ratio is sensitive to the pion sector; with pion masses that are 50% 
percent larger than the physical value, the contribution of the pion sector is drastically reduced. 
The pion sector contribution are more important for the 6th order cumulants than lower order 
cumulants. In Fig. [3] (right) we show the cumulant ratio x^lX®- The lattice results agree with 
both HRG models without pion sector or with 220 MeV pseudoscalar scalars, while the HRG 
with physical pions shows much larger discrepancy. 

3. Conclusion 

We have analyzed the fluctuations of baryon number, electric charge and strangeness in fi- 
nite temperature QCD at vanishing chemical potential. We find fluctuations and correlations of 
conserved charges are well described by an ideal, massless quark gas already for temperatures of 
about (1.5 ~ 1.7) times the transition temperature. At low temperature we find that fluctuations 
of conserved charges are well described by a hadron resonance gas up to temperatures close to 
the transition temperature. 

The current analysis has been performed with light quarks that are one tenth of the strange 
quark mass. We have shown, that higher order cumulant ratios like, for instance xf/x®' become 
quite sensitive to the pseudoscalar Goldstone mass. In numerical calculations with staggered 
fermions this also means that results become sensitive to a correct representation of the entire 
Goldstone multiplet. Calculations with smaller quark masses closer to the continuum limit will 
thus be needed in the future to correctly resolve these higher order cumulants, which will give 
deeper insight into the range of applicability of the resonance gas model at low temperature and 
the non-perturbative features of the QGP above but close to Tc. 
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